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Rapid reduced-order numerical models are being investigated as candidates to simulate the dynamics of a 
flexible launch vehicle during atmospheric ascent. There has also been the extension of these new approaches 
to include gust response. These methods are used to perform aeroelastic and gust response analyses at isolated 
Mach numbers. Such models require a method to time march through a succession of ascent Mach numbers. 
An approach is presented for interpolating reduced-order models of the unsteady aerodynamics at successive 
Mach numbers. The transonic Mach number range is considered here since launch vehicles can suffer the 
highest dynamic loads through this range. Realistic simulations of the flexible vehicle behavior as it traverses 
this Mach number range are presented. The response of the vehicle due to gusts is computed. Uncertainties 
in root mean square and maximun bending moment and crew module accelerations are presented due to 
assumed probability distributions in design parameters, ascent flight conditions, gusts. The primary focus is 
on the uncertainty introduced by modeling fidelity. It is found that an unsteady reduced order model produces 
larger excursions in the root mean square loading and accelerations than does a quasi-steady reduced order 
model. 
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= State matrix 
= Gust model scale factor 
= Generalized displacement 
= Static or mean generalized displacement 
= Generalized force due to gust 

= Generalized force due to dynamic component of gimballed engine thrust 
= Generalized force due to mean component of gimballed engine thrust 
= Fractional exponent 
= Structural stiffness matrix 
= Gust ramp width of component n, grid units 
= Structural and aerodynamic mass matrix 
= Number of lag states, Nr = 10 
= Free stream dynamic pressure, psi 
= Free stream dynamic pressure for nominal trajectory, psi 
= Generalized aerodynamic stiffness matrix 
= Generalized aerodynamic mass matrix 
= Generalized aerodynamic damping matrix 
= Unit sign function (±1) 

= Dynamic z-component of engine thrust, lbs 

= Mean ’-component of engine thrust, lbs 

= Gust time start of component n, sec 

= Free stream velocity, in /sec 

= Free stream velocity for nominal trajectory, in /sec 

= Nondimensional vertical gust velocity 

= Vertical gust reference magnitude of component n, in/ sec 
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= Nondimensional grid speed metrics 
= x location of the engine gimbal, in 
= Modified nondimensional grid speed metrics 
= Angle of attack, degrees 
= Local angle of attack, degrees 
= Local angle of sideslip, degrees 
= Roger approximation lag root matrix 
= Structural and aerodynamic damping matrix 
= Local angle of attack increment due to gust, radians 
= Matrix of eigenvectors of structural dynamics equations 

= Functional representation of structural eigenvector n projected to vehicle centerline 
= State variable matrix 
= Unsteady aerodynamic lag state 
= Structural damping ratio 
= Aerodynamic damping ratio 
= Nondimensional CFD time 
= Structural and aerodynamic stiffness matrix 


Note to Readers: The predicted performance and certain other features and characteristics of the Ares I and Ares 
I-X launch vehicles are defined by the U.S. Government to be Sensitive but Unclassified (SBU). Therefore, details have 
been removed from all plots and figures and tabulated data may have been rescaled. 

I. Introduction 

Recent interest has been directed toward developing rapid reduced-order numerical models (ROM’s) to simulate 
the behavior of flexible launch vehicles during atmospheric ascent. For example Capri, Mastroddi and Pizzicaroli 
use system identification of inviscid unsteady aerodynamics to perform aeroelastic stability analysis of the VEGA 
European small launch vehicle Silva performs system identification to extract a state representation of the unsteady 
aerodynamics of the NASA Ares I Crew Launch Vehicle (CLV) and Ares I-X Flight Test Vehicle (FTV) 2 . These 
reduced-order models linearize the unsteady aerodynamics about a mean nonlinear Computational Fluid Dynamics 
(CFD) generated flow field. There has also been the extension of this approach to include gust response, as for instance 
in the reduced-order gust model and gust response of a launch vehicle performed by Mastroddi et al. 3 . However, these 
methods are used to perform aeroelastic and gust response analyses at isolated Mach numbers. Launch vehicles 
intended for orbital insertion never remain at a given atmospheric flight condition for a significant time. What will 
also be necessary is the capability to simulate flexible vehicle dynamics through a continuous succession of ascent 
conditions. 

While the use of aircraft and launch vehicle databases at multiple flight conditions is not new (e.g. Chen et al. 4 ), it 
seems that the use of steady and unsteady aerodynamic databases to simulate launch vehicle ascent is new. The present 
paper offers one method for interpolating reduced-order models of the unsteady aerodynamics at successive Mach 
numbers. In this case, rational function approximations (RFA’s) at multiple flight conditions are used in conjunction 
with nonlinear steady databases. The interest in simulating through successive conditions is not unique. Indeed, 
there are current efforts toward combining reduced-order models at successive flight conditions, at least in the case 
of aircraft 5-7 . In the present case the requirements may be somewhat simplified by limited consideration only of 
launch vehicles. Furthermore, enhancement of methods in current industry use has certain advantages. Therefore, an 
approach is shown in which the RFA’s at key Mach numbers can be interpolated to the current flight condition. Note 
that the present simulations do not include the effect of rigid body degrees of freedom. In this paper the aeroelastic 
responses of the vehicle to gust excitation are presented. Since the engine input required for vehicle attitude control is 
one of the largest contributors to bending moment, engine thrust vectoring is used to instantaneously trim the vehicle 
at all times during the ascent. 

Slight variances in the manufacture inherently result in geometric and/or material uncertainties. If probabilities 
can be assigned, and the physics clearly understood, uncertainties in an output of interest can be estimated. It may be 
possible to assign probability distributions to design or manufacture-related quantities. Probabilistic distributions in 
structural damping can be due to varying mechanical properties of joints and material irregularities. Numerous recent 
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examples of the quantification of structural related uncertainties and risks can be cited 8-14 . Uncertainties in flow field 
quantities have also been recently investigated for a variety of flight regimes 15-18 . Probabilities derived from statistical 
atmospheric data can be assigned to flight conditions such as dynamic pressure, free stream speed or to atmospheric 
turbulence intensities and frequency. 

Uncertainties in mass and material stiffnesses can impact mode frequencies and mode shapes and ultimately the 
aeroelastic behavior of a structure. Uncertainties in aeroelastic stability, and in particular flutter onset, have recently 
been calculated for the AGARD 445.6 wing 19 . Beyond the calculation of flutter onset, probabilistic aeroelastic data 
also result in static and dynamic flight load, guidance, navigation and control and in ride quality uncertainties. Un- 
certainties can also be associated with modeling fidelity. In traditional loads analysis of flexible launch vehicles, rigid 
steady line loads are frequently used; but it is known that this method can produce results that are unconservative 20 . 
To reduce this uncertainty, various methods of higher fidelity are available, the highest being computational aeroelastic 
analysis using a CFD code. A compromise between these modeling approaches, namely a reduced-order model, is 
used in the present paper that includes the effect of unsteady flow physics. 

As an example of the use of the multi-flight condition reduced-order model, the uncertainties of several quantities 
of interest are calculated for the flexible Ares I Crew Launch Vehicle as it traverses the transonic regime. The transonic 
flight regime is chosen because it is frequently the speed range through which static and dynamic flight loads vary 
most rapidly. The aeroelastic excitation mechanism is atmospheric gusts. Here the uncertainty is considered due to a 
variety of design and flight factors that influence unsteady flexible vehicle response. The design factors are structural 
damping and frequencies. With regard to the flight conditions, the dynamic pressure, gust shapes and velocities are 
given assumed probability distributions. The parameter space encompassing all these will be simulated. The resulting 
uncertainties in the vehicle bending moment and crew module acceleration will be assessed. These uncertainties 
will be computed with an unsteady ROM and with a quasi-steady ROM, thus providing an additional assessment of 
uncertainty due to model fidelity. 

This paper will begin by discussing the implementation of gust modeling in the CFD code FUN3D. The devel- 
opment of the multi-flight condition dynamic and quasi-steady ROM’s and their verification will be discussed. Gust 
modeling for a launch vehicle ascent will briefly be discussed. Finally, these methods will be applied to uncertainty 
quantification for parameters of interest for the Ares I launch vehicle. 


II. FUN3D Code 


A. Basic Solver 

The Navier-Stokes code used in this study is the Fully Unstructured Navier-Stokes Three-Dimensional FUN3D code, 
a finite-volume unstructured CFD code for either compressible or incompressible flows 21 - 22 . Flow variables are stored 
at the vertices of the grid. FUN3D can solve the discrete compressible Euler or Reynolds-averaged Navier-Stokes 
(RANS) flow equations either tightly or loosely coupled with a turbulence model on mixed element grids, including 
tetrahedra, prisms, pyramids and hexahedra. In the present study the RANS solver and the loosely coupled Spalart- 
Allmaras turbulence model 23 are used on an all tetrahedron grid. FUN3D employs an implicit upwind algorithm in 
which the inviscid fluxes are obtained with a flux-splitting scheme. The low dissipation flux splitting scheme for 
the inviscid flux construction, and the blended Van Leer flux limiter were used 24 . The full Navier-Stokes viscous 
terms are used 21 . Steady state solutions are accelerated to convergence by the use of local time stepping 25 . Domain 
decomposition enables distributed computing on multiple processors. 


B. Field Velocity Method of Modeling a Gust 

Recent innovations have allowed the simulation of a gust with a CFD code using the Field Velocity Method (FVM) 8,26-28 
The FVM physically introduces gusts into a CFD code by utilizing grid velocity. 28 Normally grid velocity would be 
associated with physically moving the grid. However, it is possible to introduce an arbitrary perturbation velocity in 
a stationary grid by prescribing the grid velocity at either all or some of the field grid points without actually moving 
the grid. If the gust shape is defined by a streamwise variation in a perturbation velocity, expressed as 

w g (x,t) =f(t — x/Uco) , t>x/Uoc (1) 
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the gust profile would be introduced into the flow field by modified grid speed metrics 

Xxi + yxj J rZzk = x x i + yxj +{zT — w g )k . ( 2 ) 

To date this method has been used for relatively simple to moderately complex configurations. 29 Further use and 
development of the method have been performed by Singh and Baeder 8 , Zaide and Raveh 27,30 , Raveh 29-31 and Wang 
et al. 32 . 

This approach has been implemented in the FUN3D CFD code. Gust profiles such as sharp edged or compound 
sequences of one-minus-cosine increments are used in the present paper. A CFD calculation incorporating a gust in 
this fashion will be compared with a gust introduced into the dynamic reduced-order launch vehicle model. 


III. Reduced Order Models 


In order to clarify what is new, the discussion here will focus on the flight condition interpolation, gust input and 
response, and engine trim. The details of the method for isolated flight conditions are discussed elsewhere 20,33,34 . In 
reference [20], the solution is partitioned into linearized dynamic and nonlinear quasi-steady mean components. The 
dynamic solution is linearized about the nonlinear mean aeroelastic solution. The mean solution is obtained using 
equation 1 1 of reference [33] and steady state CFD generated aerodynamic line loads at Mach numbers 1 .00, 1 .20 and 
1.55. Line loads are linearly interpolated to ascent Mach number, a / and /!/ to compute the generalized forcing. The 
quasi-static mean equations can be written 

[/] g s " + [2co\ gs + [or] fe} = {G, } + {Ges} (3) 


where G s is the generalized forcing due to aerodynamic line loads. The primes in equation 3 refer to differentiation 
with respect to a pseudo time used to converge the static mean problem. This paper also includes generlized force G es 
due to the mean transverse gimballed engine forcing, where for mode n, G esn = T- v <t>„ (x g ) . Engine thrust is used to 
trim the mean component of the vehicle pitch moment. 

The dynamic response of the vehicle is computed by linearizing the lineloads around the local mean a/, /3/ and 
static generalized force G s . A state space model of the aeroelastic system with lag states (<*) is developed by using 
RFA theory. The state space model is given by 


where 


and 


{i} = [A]{x} + {G g } + {G ed } , {x}={g 1 g^} T 
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[M(M 00 )\=I-p 00 [Q m {M 00 )\ , l£l(M„)]=[(0 2 ]-p oc Vj[Q k (M oo )\ , [A(Hf»)] = [2(>] -pcV* (M„)] . 

(6) 

G ed is the component of engine thrust required to trim the gust forcing, for mode n given by G edn = T zd < t>„(x g ). G g 
on the right hand side of equation 4 is a vector of generalized forces due to gust. Gust modeling will be discussed 
in the next section. Through the aerodynamic mass ( Q m {Moa)), damping and stiffness ( Qk(Moo )) matrices 

and dynamic pressure the matrices M(Mj), il(M r „) and A (AT,.) are functions of Mach number. The line loads 
also vary with Mach number. Because of propellant burn, structural frequencies (to) and mode shapes are functions of 
time into ascent. As discussed in reference [20] modes 1 and 2 are modeled in a fully unsteady sense using an RFA 
approximation. The unsteady aerodynamic submatrix is obtained via the Roger approximation 33 . This approximation 
results in equations 12 and 13 of references [20,36]. The remaining modes are modeled with quasi-steady aerodynam- 
ics using equations 14 and 15 of reference [33]. Definitions of terms in these equations are found elsewhere 33 . In this 
paper, the structural and aerodynamic equations at Mach numbers 1.00, 1.20 and 1.55 are used to simulate the ascent 
of the vehicle through this Mach number range. Based on trajectory data, the steady and unsteady databases at these 
Mach numbers are interpolated to each time step. The constant matrix e is defined in references [20,36]. 

To create a system that can traverse several Mach numbers in a stable manner, it was found essential to enforce 
the constraint that the same values of y are used for all Mach numbers. Therefore, 10 lag roots at Mach 1.00 were 
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obtained and the same used at all Mach numbers, 1.00, 1.20, and 1.55. Using the same Mach 1.00 lag roots ( 7 ), 
unsteady aerodynamic models ( d , Qk, Q%, and Q m ) were obtained at Mach numbers 1.20 and 1.55. With the complete 
database assembled, the numerical integration of equation 4 was performed with a 5 ,h order implicit Euler backward 
difference scheme. The line loads, mode shapes, modal frequencies and unsteady aerodynamic data at each time step 
are linearly interpolated between the data at the three Mach numbers. 

The response is also calculated using linearized quasi-steady line loads. Linearization around the local mean a/, 
Pi results in a state space model given by 

ix} = [A]{x} + {G g } + {G ed } , { X } = {g,g} T (7) 

where 


The definitions in equation 6 still apply, but the aerodynamic damping and stiffness matrices are given by equations 
14 and 15 of reference [33], and the aerodynamic mass matrix is zero. 

Because separate MSC.Nastran™ eigenvalue analyses were performed on the structural models at each of the 
three flight conditions, the mode shape vectors were assigned differing signs by the software. Since it is possible to 
define 

[,n] = mm T , [K) = m[co 2 ][p} T (9) 

the sign on eigenvectors is arbitrary. Accordingly, the signs of the mode shapes at Mach numbers 1.20 and 1.55 were 
changed so that mode 1 maintained the same sign orientation at all Mach numbers. This sign change too was necessary 
for a stable multi-Mach number simulation. Finally, what is expected to be a conservative approach to engine loading, 
engine gimballing is used at each time step to instantaneously trim the vehicle response to gust loading 37 . 


0 I 

-AT 1 A 


( 8 ) 


IV. Gust Model 

In this study composite one-minus-cosine gust profiles typical of the Statistical Discrete Gust (SDG) method will 
be used. This method assumes that gusts can be built up as an assembly of discrete, smoothly varying ramp-hold 
functions 38 - 39 . Four one-minus-cosine ramps are used separated by hold functions. In the present analysis, gusts are 
modeled by adding a local angle of attack increment due to gust, given by 

1 ’ 

~ 52 w ref n {Lg n /b g ) Sgn„(l — cos(G„))/Uoo (10) 

z M=1 

where 


n(t-t gn - x /U aa ) 

for t gn + x/Uoo < t < t gn +X/U 00 + L gn / U 0 0 

( 11 ) 

Lg n /Uoo 

II 

for t >t gn +x/Uoo+L gn /Uoo ■ 

( 12 ) 


Figure 1 shows a typical velocity profile of a composite gust with 4 segments and the phase and width parameters 
t g and L g for segment 3. The resulting local angle of attack increment due to gust produces a loading increment, 
computed from the steady rigid line loads database. Integration over the vehicle length produces the generalized force 
increment G g due to gust. This provides the forcing on the right hand side of equations 4 and 7. The response to the 
gust is found in the evolution of the generalized variables and the lag states. 


V. Verification of Performance 

Figures 2-4 show the aerodynamic dampings and frequencies resulting from the present unsteady ROM and the 
FUN3D computational aeroelastic code separately computed at each of Mach numbers 1.00, 1.20 and 1.55. Note 
that aerodynamic damping as defined here is due to the motion of the aeroelastic modes and is a function of dynamic 
pressure. The unsteady ROM solutions are computed at q nom and \32q nom while time marching FUN3D solutions are 
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shown at q nom . Most of the solutions are at a = 0.0 degrees except the Mach 1.55 FUN3D time marching solution, 
which is at a = 2.5 degrees. (The aerodynamic damping at Mach 1.55, a = 0.0 and a = 2.5 degrees are similar.) 
In each case the frequency and aerodynamic damping values in these figures were obtained from a log decrement 
of the time response to an initial generalized modal velocity. Note that at Mach numbers 1.00 and 1.20 mode 1 
is aerodynamic ally undamped and becomes further undamped at higher dynamic pressure. At Mach 1.55, mode 1 
aerodynamics is damped and becomes further damped at higher dynamic pressure. With the exception of Mach 1 .20 
where the unsteady ROM model slightly over predicts mode 2 damping and mode 1 undamping, the unsteady ROM 
results compare very well with the time marching computational aeroelastic simulations. This initial verification 
indicates that individual ROM databases are consistent with the full order CFD solutions. The progression from 
negative aerodynamic damping in mode 1 to positive damping at higher Mach numbers is a feature that will be 
important in the simulations to follow. 

Since a gust model is included in the reduced-order analyses, it is of interest to compare the gust responses of the 
reduced-order models with that from the FUN3D time marching gust simulation. Figures 5-6 show the generalized 
displacements and force time histories comparing Ares I CLV gust simulations using the full order FUN3D CFD and 
unsteady ROM methods. Figure 5 shows the Ares I generalized displacement and force time histories due to a sharp 
edge gust equivalent to a 2 degree angle of attack increment. In this case critical structural damping is used to damp out 
structural oscillations. Because of the time lag required for lift to be generated, there is a period of time over which the 
steady state displacement and forces develop. Although there are minor differences, the full order and reduced-order 
solutions compare quite well. 

Figure 6 presents a somewhat more challenging case. The composite gust profile with four one-minus-cosine 
gust ramps, and associated time sequencing produces the response shown. The reduced and full order generalized 
displacement time histories match very well. The generalized force requires somewhat more discussion. The quasi- 
steady and unsteady ROM’s include the aerodynamic modal forcing implicitly in matrix [A] of equations 5 and 8. The 
gust forcing is included explicitly in the term { G„ } . Therefore, the unsteady ROM generalized force shown in each 
of Figures 5 and 6 are only the gust forcing (G„ in equations 4 and 7) while the full order FUN3D generalized force 
is due to all aerodynamic effects combined. Additional dynamics are therefore evident in the FUN3D time accurate 
solution that are not seen in the gust forcing increment in the ROM solution. But the displacements and forces give at 
least a qualitative agreement. 

Having generally verified the performance of the model for individual flight conditions and the gust simulation 
methodology, the progression through flight conditions is now demonstrated. The mode 1 time history of the response 
to an initial generalized velocity excitation is presented in Figure 7. Here, a comparison is shown between the unsteady 
and the quasi-steady ROM’s. These simulations traverse the entire Mach number range from 1.00 to 1.55. Note that 
in the unsteady ROM the mode is initially undamped. Beyond Mach 1 .20 mode 1 in the unsteady ROM simulation 
becomes aerodynamically damped. Comparing these results with the dampings of Figures 2-4, the aerodynamic 
damping of mode 1 at each Mach number is correctly reproduced in the unsteady ROM simulation. Also shown in 
Figure 7 is a quasi-steady simulation over the same Mach number range. The quasi-steady simulation results in 
a strongly damped mode 1. This result is of course physically inaccurate, in that the quasi-steady approach does 
not include the dynamics resulting from the unsteady aerodynamics and associated time lags. The result is, though, 
consistent with the quasi-steady results shown previously 34 . 

Finally the dynamic thrust gimballing used here is demonstrated in Figure 8. Note that as discussed earlier, the 
mean part of the solution is trimmed for changes in Mach number and dynamic pressure in a static sense. On the other 
hand, the linearized dynamic solution introduces the gust excitation into the solution. The gust excitation produces an 
instantaneous dynamic thrust angle change to trim at each time step. An instantaneous thrust angle change to trim is 
physically unrealistic since a control system will not fully trim the vehicle at each point in time. The present method 
however, can be considered conservative in that it will produce a higher bending moment than dynamic thrust control 
would produce 37 . Figure 8 shows the gust velocity excitation and the dynamic thrust gimbal angle resulting from that. 
The bending moment time history is also shown. The bending moment shown is measured at the streamwise midpoint 
of the vehicle. The introduction of thrust gimballing is clearly seen to increase the peak dynamic bending moment. 


6 of 16 


American Institute of Aeronautics and Astronautics 



VI. Example - Quantification of Uncertainty 


The quantities of interest are the uncertainty in the maximum bending moment, rms bending moment at the stream- 
wise midpoint of the Ares I launch vehicle, and rms acceleration at the Crew Exploration Vehicle (CEV). To compute 
uncertainties in bending moment due to variations in flight and design parameters, a total of 32,000 computations have 
been performed using the quasi-steady and 32,000 using the unsteady ROM’s. These solutions comprise all possible 
combinations of the chosen flight, design and gust parameters. Each simulation was exposed to a gust profile com- 
posed of four one-minus-cosine shapes spaced out by four time start parameter values. Three sets of gust shape signs 
were used. The gust time starts (t gn ,n = 1,4) were varied between three sets of four values. For maximum launch 
vehicle loading due to gust, a 350 in/sec one-minus-cosine wave form was tuned to the lower frequency modes by 
setting the maximum value of L s„- Standard atmospheric data suggests this will give a 99 percent closure and 90 
percent confidence level loading 40 . The gust magnitudes (w re f n ,n = 1,4) take on three sets of four values, between 
100 and 350 in/sec. The gust width parameters (L g ,n =1,4) took eight sets of four values variously spaced between 
minimum and maximum one-minus-cosine lengths. 

Computations with all of these parametric variations were performed both for the purely quasi-steady and the 
unsteady ROM’s. Furthermore, to assess the likelihood of a given maximum bending moment value, Gaussian prob- 
ability distributions centered about a nominal value have been assigned to each of the parameters. Gust velocities 
are also assigned probabilities reflecting the atmospheric data of reference [40], (For instance, a 350 in/sec gust is 
given here a probability of approximately 1 percent.) In this paper the total probability of an event is the product of 
probabilities of each of the design, flight and gust parameters. 

Illustrative of the results obtained by these simulations. Figure 9 shows the total bending moment time history for 
the most extreme condition, computed with both the quasi-steady and the unsteady ROM’s. This condition is the case 
giving the largest maximum total bending moment. Total bending moment is the combination of the static mean and 
dynamic bending moment components. The dynamic component is caused by the gust excitation plus required thrust 
gimballing to trim during gust. The parameters resulting in this case are the highest dynamic pressure, lowest first 
mode frequency, lowest structural damping, highest gust velocities with best tuned width. Note that for this case, the 
dynamic component dwarfs the static mean component of bending moment. The gust excitation is completed within 
the first 2 seconds of the simulation while the flexible response continues much longer. The largest absolute magnitude 
of bending moment in the quasi-steady simulation occurs at about 2.5 seconds while the largest absolute magnitude 
of bending moment in the unsteady simulation occurs at about 9 seconds into the simulation. Because of the unsteady 
aerodynamics of mode 1, the unsteady ROM results, over the entire course of the simulation, in larger excursions 
of bending moment than the quasi-steady simulation. The response of mode 1 with the quasi-steady ROM remains 
damped through out. Even though the unsteady ROM gives a 15 percent larger maximum bending moment than does 
the quasi-steady ROM, the probability of the combination of parameters producing this result is extremely small, well 
beyond a 99 percent closure requirement. 

Figure 10 presents outlines of the histograms of the maximum total bending moment for every combination of 
design, flight and gust parameters. The resulting values are assigned among 150 bins. Typical histograms show 
bar plots, however, in the present case, bars would obscure the comparison of the two sets of results. Therefore, 
only the outlines of the histograms are shown. Summed probabilities for all occurances of each maximum bending 
moment value form the outlines of the histograms shown. The quasi-steady and unsteady ROM’s have nearly the same 
distribution of maximum bending moment. As just discussed, only in the most extreme and most improbable cases is 
there a significant difference in maximum bending values. The 95 percent values of maximum bending moment are 
nearly identical for the two methods. The 95 percent value of the rms bending moment, as indicated in the histogram 
outlines of Figure 1 1 computed with the unsteady ROM is significantly higher than that computed with the quasi- 
steady method. This is indicative of the greater persistance of the mode 1 oscillation for the unsteady ROM. Figure 12 
shows the histogram outline of the rms ^-acceleration at the CEV. The persistance of mode 1 dynamics is also evident 
in the higher 95 percent acceleration level experienced at the crew module. 

VII. Concluding Remarks 

This paper has presented an approach to simulating flexible launch vehicle ascent through multiple flight condi- 
tions. The approach includes a simultaneous computation of a nonlinear mean solution as well as a linearized dynamic 
solution. Although currently restricted only to flexible vehicle modes, the approach is applicable to rigid body degrees 
of freedom as well. The approach shown here is verified by comparing aerodynamic damping at three transonic Mach 
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numbers. A gust model is also included in the reduced-order model, and verified by comparison with a direct CFD 
simulation of a gust. Dynamic and quasi-steady reduced-order model databases are constructed for three Mach num- 
bers through the transonic regime. The dynamics of the flexible vehicle is compared for the two methods. An example 
of the use of these reduced-order models and multi-flight condition databases, the uncertainty in the maximum bending 
moment, rms bending moment at the streamwise midpoint of the Ares I launch vehicle, and rms acceleration at the 
crew module are computed for ascent traversing the transonic regime. The resulting probability distributions show 
that the unsteady reduced-order model results in higher 95 percent rms vehicle bending moment and crew module 
acceleration than the quasi-steady reduced-order model. 
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(a) Unsteady ROM. 


(b) FUN3D time marching, q„ om . 


Figure 2. Aerodynamic damping versus frequency, Mach 1 .00, a = 0.0 deg.. 
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Freq ~ Hz. 



^a C 


(a) Unsteady ROM. (b) FUN3D time marching, q nom - 

Figure 3. Aerodynamic damping versus frequency, Mach 1.20, a — 0.0 deg.. 


(a) Unsteady ROM, a — 0.0 deg.. (b) FUN3D time marching, a = 2.5 deg., q nom - 

Figure 4. Aerodynamic damping versus frequency, Mach 1.55. 
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FUN3D time accurate gust 
unsteady ROM 




Figure 5. Generalized displacement and force time histories due to a sharp edged 2 degree gust, Mach 1.00, a = 0.0 deg.. 




Figure 6. Generalized displacement and force time histories due to composite one-minus-cosine gust, Mach 1 .00, a = 0.0 deg.. 
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Figure 7. Mode 1 generalized displacement time history due to initial generalized velocity. 
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(a) Gust velocity. 


(b) Dynamic component of engine gimbal angle. 
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(c) Dynamic component of bending moment. 

Figure 8. Effect of thrust gimbaling to trim. 
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Figure 9. Time history of total bending moment. Most extreme loading case. 



Max Bending Moment 


Figure 10. Probability distribution of maximum total bending moment. 
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Figure 11. Probability distribution of rms total bending moment. 



( a M ms 

Figure 12. Probability distribution of rms z-dir. acceleration at CEV. 
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